Thermodynamic functions within the van Hove BCS model: Symmetry mixing effects. 
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All the new layer perovskite superconductors seem to show a phenomenon of symmetry mixing 
with repect to the order parameter. An analysis of the different alternative of mixing and how far 
could them be presented is carried out. For the particular case of s+id symmetry of the gap, the 
temperature dependence of the specific heat (C e3 ) and the thermodynamic critical magnetic field 
(H c ) are calculated. A double peak transition is observed on C es (T) in the mixed regime while the 
single peak behavior is recovered for a purely symmetric state ( s or d). C es presents a quadratic law 
at low tempeartures for a d-wave gap and for the s-wave one the typical exponential attenuation. 
The temperature dependence of H c shows a clear phase transition of second order at temperatures 
where the d-wave component becomes negligible. A comparison with other results and experiments 
is done. 
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I. INTRODUCTION. 

<N 

The discovery of the high temperature superconductivity (HTSC) in layered perovskites has stimulated many 
, theoretical studies that range from phenomenological approaches to microscopic models. Some of them follow the 
traditional Fermi liquid picture with a BCS |ffl] or Eliashberg formulations, while others, often focusing on the 
, new electronic fluids, use to follow the so called marginal and nearly antiferromagnetic Fermi liquids [|J. However, 
up to date a fully consistent interpretation of the large amount of existing experimental data is not found in the 
literature and further efforts are needed for the understanding of the pairing mechanism in cuprate perovskites. In 
this context, it is important to compare the thermodynamic functions of these new compounds and the old layered 
superconductors with low critical temperatures, looking for similarities and differences between them. 

The cuprate perovskites are layered materials with highly anisotropic properties which show an almost dispersionless 
electronic bands with flat sections in the vicinity of the Fermi level (sf). This leads to the appearance of the so called 
van Hove singularities (vHs) in the electronic density of states (EDOS) Q. On the other hand, in these compounds 
the possibility of gap symmetries different from the s-wave one is a matter of debate up to now MM. One of the 
first attempts to explain the great variety of behaviors experimentally observed was the BCS model within the 
vHs scenario using the s-wave symmetry for the order parameter j8| |Ti[. Nevertheless, several experiments, such 
as tricrystal corner junction tunnelling in Y Ba2Cu^Oj-s |l2| show strong evidence that the gap symmetry was 
dominated by a d x i^ y i component. Following these ideas and in the framework of the BCS model, in Refs. |]l3| , |I^| , 
the temperature dependence of the gap, the ratio 2A(0)/&bT' c , and the specific heat C es as a function of temperature 
were calculated. While this approach seems to yield results more close to the observed experimental behavior of those 
magnitudes, several reports ( see j(| and references therein) are very hard to explain in terms of purely symmetry 
states. In this sense, recent experiments on c axis tunneling in the Pb — YBCO system | fl5| (and references therein) 
reveals that the superconducting order parameter should have a small s wave contribution in addition to the d wave 
one. With that aim, in Refs. [[l5|-|i"8| the phase diagram for the s + d and s + id gaps as a function of the filling factor 
and the driving parameter Vd/V s were studied in the framework of the BCS model. The purpose of this paper is to 
study the effects of the symmetry mixing of the gap on the thermodynamic functions as the electronic specific heat 
C es and the critical magnetic field H c within the weak coupling van Hove BCS (v-BCS) model. 

The paper is organized as follows. In Sec. 2 the theoretical background and gap symmetries are introduced within 
the van Hove singularity scenario. In Sec. 3 the thermodynamic functions of the superconducting material as a 
function of temperature are obtained. Finally, in Sec. 4 the main conclusions of the paper will be given. 

II. TOY MODEL FOR A MIXTURE OF GAP SYMMETRIES. 

In agreement with the layered structure of the new superconducting materials we start from a tight binding for- 
mulation of the band structure with the following dispersion law: ej: = —2t(cosk x a + cosk y a), where t is the coupling 



f 



parameter between the Cu and O atoms, k x and k y are the components of the wavevector, and a is the lattice constant. 
This band model leads to a density of states similar to that obtained in Ref. Furthermore, we have considered 
the weak coupling limit of the BCS theory including gap symmetries different from the often used s- or d- waves as 
well as its possible mixing. We are discussing mainly the mixed s + d and s + id cases which seem to be closest to 
experimental results. The tetragonal lattice symmetry assumed in the band structure will impose restrictions on the 
possible forms of and for the gap the following expression are assumed [HfjJ : 

Ag(T) = A|(T) + CTp (id)A|(T), (1) 

with 



A|(T) = A s (0) ff (T/T c s ), (2) 
Ag(T) = A d (0) (cos(k x a) - cos(k y a))g(T/T?). (3) 

The term exp(i#) in Eq.([j]) enables one to fix the relative phase between the s and d symmetries. For 9 — we 
are dealing with the s + d mixing for which the order parameter is a real quantity, while for 9 = ir/2 the s + id 
mixing a complex order parameter is presented. On the other hand, two different critical temperatures, Tjf and 
should appear in agreement with the respective symmetries of the gap, so that the critical temperature of the whole 
superconductor will be the maximum among these two values (T c = max(T c s , T^)). Which of the two symmetries is 
responsible for the critical temperature depends on many factors p Ul^JT^I as for example, at higher filling values the 
d x 2_ y 2 gap is favored, while for lower filling factor we have a pure s-gap Jig] . Irrespective of the conditions favoring 
one or another gap symmetry mixing we are interested on the influence of this on the thermodynamic functions. With 
this purpose and to simulate how the material evolves from s to d-wave, the following equations for T c are used: 



k B T c s (c) = 1 

2t exp[(c-0.3)/0.04] + l 
k B T c d (c) 1 1 



(4) 



2t exp[(c- 0.5)/0.06] + 1' v ; 

where c is a parameter that drives the amount of d-wave component and a suitable estimation for this parameter 
could be the ratio Vd/V s . The values c = and c = 1 correspond to the nearly s- and d-wave gaps, respectively. 
This model is consistent with results obtained in Refs. |l6| . |l8| where three regions are well defined. Two of them have 
a dominant gap symmetry, s or d and the third one a mixture of these two symmetries. On the other hand, and 
irrespective of the gap symmetry assumed, we denote by <?(T/T^) the temperature dependence of the gap which with 
a good accuracy can be represented by the following relation |24| : 



T 

9{jt) = tanh 



T c -T 
T 



(6) 



It is worth noticing that the function g(T/T c ) differs slightly from the correct BCS dependence and the slope at T c is 
about 1.6 times the true one. 

Since the thermodynamics of any superconductor is governed by the concrete topology of the quasiparticle excitation 
energy, in Fig.l the contour lines generated by the surface E^ — ^/e| + | A^| 2 are plotted for the s-wave (upper left) 

and d-wave (upper right) pure symmetry cases as well as the s + d-wave (bottom left) and s + id-wave (bottom right) 
mixing at zero temperature. The contour lines of the symmetry mixing cases are plotted considering the d-wave 
component as the dominant one and assuming T c s = 0.7T c d . In all cases, E^ is expressed in units of 2t as a function of 
the dimcnsionless wave vectors components X — k x a and Y = k y a. As can be seen from these pictures, the function 
does not have zeros in the pure s-wave case while for the d-symmetry isolated zeros at (±7r/2, ±7r/2) are present. 
These zeros lead to qualitatively different results in the behavior of thermodynamic quantities at low temperatures. 
While for a s-wave gap an exponential attenuation in C es is expected in agreement with the activation character of 
the quasiparticle spectrum, for the d-wave case a power law follows, due to the absence of any threshold. Concerning 
the mixing of symmetries, it is clear that the s + d and s + id gaps should lead to a different character form of the 
thermodynamic functions. Since the symmetry of the contours associated with the s + d mixing is lower than the 
tetragonal one this kind of mixture is more suitable for lattice with a small ortorhombic distortion, as for example in 
the YBaiCu-sOi-y system. 
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III. THERMODYNAMIC FUNCTIONS. 



From experimental point of view a double peak transition in the specific heat (C es )has been observed for YBCO 
Jl9| plj. In Rcf. Q a double step at 88 K and 92 K were reported by adiabatic continuum method. Moreover, in 
Refs. |l9|,pT|] a review of the experimental evidences about the double transition in C es for ReBa-iCuzOi (Re = Y, Gd) 
system has been presented. These experimental facts are in correspondence with the mixed gap scenario. Starting 
from the above motivation in this section we are dealing with the specific heat in the superconducting state and the 
thermodynamic critical magnetic field taking into account the gap symmetry mixing. 

The electronic specific heat in the superconducting state can be written as: 

k 

f(EjA being the Fermi-Dirac distribution function. The derivative of the square of the gap in the case of a mixture 
of symmetries can be cast into: 

d|A r (r)| 2 dAUT) 2 dAi(T) 2 dA'JT) . dAf(T) 

Following Eq.(Q) and irrespective of the gap symmetry, the derivative of A(T) and A 2 (T) with respect to the tem- 
perature are given by: 



; if T = T C , 



2T r I — f — A2 (°) tanhpyS^^ 1 



2A 2 (0) 

Tr 



if T = T r 



As can be seen from Eqs.(||) and (||), in the case of s + d mixing an square root singularity in the specific heat is 
obtained if Xj? and are different. Otherwise, a smooth dependence on T/T c is achieved. However, up to now 
experimental reports on the presence of this kind of singularity are still absent. Hence, in order to observe the s + d 
symmetry the critical temperatures for the both components (A s and A^) must be the same. It is important to 
remark that this result have nothing to do with the EDOS or any criteria of stability for the phase under certain 
crystalline symmetry. For the s + id mixing the square root singularity does not appear, because the relative phase 
between the two components of the gap is equal to ir/2 and the cosine function in Eq.(||) is zero. Further we shall 
deal with the influence of the s + id mixing on C es (T) and H C (T). 

Transforming the sum into integrals over the first Brillouin zone and applying similar methods as in Ref . Jig] , C es 
can be cast into: 

c „ ( r,= « (»,>&> r « r*»("0(i-/W) 



where 



VT - 7? k B T c ' K T' J_ Uc J V[ W(v,Q 
^.0-^)^[(^r) a + (^rW«,0]], (H) 



2tt 2 ^(1 - v)(l + v)(v - vl)(v2 - v) 

D(v,0 = - (d + 2v) vl = -(l + 0, u2=(l-0, 

R is the universal gas constant and io c — e c /2t is the dimensionless energy cut-off. In the following for 2A S (0) /fesX^ 
and 2Ad(0)/kBT c 1 the values of 4 Q and 4.26 |13] have been used, respectively. The integration limits over v (lower 
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vi and upper v u ) are defined as follows: for — lu c < £ < 0, = vl and v u = 1, while for < £ < w c , vi — — 1 and 
v u = v2. For numerical evaluation, the values ui c = 0.07 (e c = ZhmeV) and ksT c /2t — 0.017 (T c = lOOif) have been 
used. The function N(£,v) is a measure of the EDOS for the respective symmetries of the gap as it can be seen by 
the following integrals ]is|] : 

N d (0= f ' dvN(v,£) D 2 {v, £), (12) 

J Vi 

N s (0= dvN(v,0- (13) 
N s (^) coincides with the exact EDOS of the tight binding model on the square lattice, given by a complete elliptic 



integral of the first kind (see Ref. |11|) and Nd(C) could be regarded as an effective EDOS for the d-wave pairing. 

The dependence of C es /C en (C en is the specific heat in the normal state) on the reduced temperature T/T c is 
presented in Fig. 2 for c = and c = 0.96. The first case corresponds to a nearly pure s-wave while the second has a 
d-wave character. As can be seen from this picture the jump due to the s-wave gap at T c is smaller than that of the 
d-wave one. This is in contrast with the conventional anisotropic superconductors, where the specific heat jump in 
the isotropic case is larger than that for the anisotropic one p3]. The above fact is explained as follows. Taking into 
account Eqs.(||), (||), and (||) the specific heat jump can be written as: 



AC(T C ) 



18.14fc fl f_ 2 di N d (0 [1 - /(Id)] ; d-wave 

16k B f_ 2 d£ N s (0 /(Id) [1 - /(|d)] ; s-wave 



We note that the term / (1 — /) is the same for both symmetries. The square of the cosine difference appearing in 
the EDOS Nd(£,) for the d-wave case plays a key role for the understanding of the difference of AC(T C ) between the 
two symmetries. It is important to note that -/V^(£) function has a stronger singularity than that coming from the 
s-wave gap explaining why the AC(T C ) value for d- waves is larger than the corresponding to the s-symmetry. In the 
inset a comparison between the s- and d-wave specific heats at low temperature is shown. 

The evolution in term of the driven parameter c for the temperature dependence of C es (T) /C en (T c ) is shown in 
Fig. 3. For c < 0.385 the s-wave gap is dominant (see Fig. 3(a)), while for c > 0.385 the T c is ruled out by the d-wave 
(see Fig. 3(b)). The double peak transition is better seen for 0.37 < c < 0.40. In Fig. 3(a) (Fig. 3(b)) the peaks 
numbered by 1 to 4 (1 to 3) correspond to the d-wave (s-wave) contribution. The results depicted in Figs. 3(a) and 
(b) constitute a direct consequence of the difference in the jump AC(T C ) for the mixed symmetries. These results 
are qualitatively in good agreement with those reported in |pj|-^l|| where a double peak transition has been observed. 
It is worth noticing that in recent measurements on Y , s Ca ,2Ba 2 Cu 3 Oj-s [§3| this structure apparently seems to 
be absent. Nevertheless, these data show, for certain values of Oxygen concentration, a small feature that could be 
interpreted in the framework of the mixed waves scenario. 

The behavior of C es (T)/C en (T c ) at very low temperatures is depicted in Fig. 4 for values of c > 0.4, where the 
d-wave contribution is dominant. The quadratic behavior observed in a wide range of temperatures for the d-wave 
case is in contrast with the linear dependence obtained in Ref. |Q. This quadratic law is consistent with the 
experimental measurements of Ref. p5j and the phenomenological approach which assumes a linear dependence on 
energy for the EDOS in the superconducting state. As soon as c reaches the values 0.52, in the low temperature 
region C es (T)/C en {T c ) goes through a small and broad shoulder around T = 0.005T C superposed onto the quadratic 
behavior law. This anomaly seems to be related with the concrete topology of the gap directly reflected in the lines of 
constant quasiparticle energy of the d-wave superconductor. It appears only when the van Hove and the d-wave 
gap are simultaneously present, otherwise they can not appear. 

The thermodynamic critical magnetic field as a function of temperature can be obtained from the usual BCS 
expression Q. In Fig. 5 the dependence of the reduced critical magnetic field H C (T) / H c (0) on the reduced temperature 
T/T c is shown for different values of c parameter. As can be seen see in Fig. 5(a), for c < 0.39 the s-wave component 
is dominant and the phase transition to a single symmetry state is clearly observed. Here, the slope at T c associated 
with the s-wave gap is lower than that coming from the conventional BCS. Moreover, when the d-wave component is 
dominant the phase transition is almost unobservable and the temperature dependence in this case shows an small 
departure from the old parabolic law (compare Figs. 5(a) and (b)). This quasi-parabolic law in the d-wave case is 
a direct consequence of the zeros in the quasiparticle excitation energy, which also provides a power law in the low 
temperature regime of the specific heat (see inset in Fig. 2). 
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IV. CONCLUSIONS. 



A toy model within a van Hove BCS scenario which takes into account the existence of mixing of gap symmetries 
has been developed. Within this framework the specific heat dependence on the temperature was computed resulting 
in a double peak transition whose resolution is very sensitive to the amount of the d-wave component present in the 
gap. For a pure d-wave gap, a quadratic law at low temperatures with a superposed shoulder at T = 0.005T C is 
obtained. On the other hand, AC jump at T = T c of the pure d-wave component is larger than that coming from an 
s-wave one. These features are due to the joint effect of the gap anisotropy and van Hove singularity. It is expected 
that the influence of anisotropy will be reinforced in those directions where it is allocated. Our calculations show 
that the critical magnetic field follows a quasi-parabolic law with temperature which results from the isolated zeros in 
the quasiparticle excitation energy. The anisotropy effect along with the van Hove singularity establish a qualitative 
difference with respect to old conventional superconductors where a constant EDOS is considered. Moreover, it seems 
that a model with a dominant d-wave contribution plus an small s-wave component may show a closes fitting with 
the reported experimental measurements. 
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Figure captions 



Fig.l Lines of constant quasiparticle excitation energy in unit of 2t as a function of the dimensionless wave vector 
components X — k x a and Y = k y a for a gap with a pure s-wave (upper left corner), and d-wavc (upper right 
corner) symmetries, an s + d (lower left corner) and s + id (lower right corner) mixing. 

Fig. 2 Electronic specific heat in the superconducting state in units of C en (T c ) as a function of the reduced temperature 
T/T c for the values of c = (s-wave) and c = 0.96 (d-wave gap). The inset shows the range of low temperatures. 

Fig. 3 Electronic specific heat C es (T)/C en (T c ) as a function of T/T c for different values of the driven parameter c in 
the region of gap symmetry mixing: (a) c = 0.37 (1), 0.375 (2), 0.38 (3), and 0.385 (4); (b) c = 0.39 (1), 0.395 
(2), and 0.4 (3). 

Fig. 4 Very low temperature behavior of the superconducting specific heat for different values of c: 0.43 (1), 0.46 (2), 
0.49 (3), 0.52 (4), 0.55 (5). A superposed shoulder is clear at T = 0.005T C for values of c equal to 0.52 and 0.55. 

Fig. 5 Thermodynamic critical magnetic field in units of H c (0) as a function of the reduced temperature T/T c for 
different values of the driven parameter c in the region of gap symmetry mixing: (a) c = 0.37 (1), 0.375 (2), 
0.38 (3), and 0.385 (4); (b) c = 0.39 (1), 0.395 (2), and 0.4 (3) (see text for explanation). 
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